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Recently, Onari and Kontani submitted a paper 1 which 
criticizes our recent theoretical study 2 - on the neutron 
scattering experiment as a probe for determining the su- 
perconducting gap in the iron pnictides. In their paper, 
Onari and Kontani have developed a formalism in which 
the imaginary part of the dynamical spin susceptibility 
(Imx s (q, u))) in the superconducting state can be more 
accurately calculated especially in the oj < 2A regime, 
where A is the superconducting gap. In section IIIC of 
their paper, they mention that the conclusions of our 
paper are "incorrect based on inaccurate numerical cal- 
culation". In the present Comment, we show that this in 
fact is not correct. 

First of all, Onari and Kontani emphasize the differ- 
ence of the calculation method used, but actually, a more 
important difference between their calculation and ours 
lies in the parameter values taken. The most important 
difference is while we take the strength of the quasipar- 
ticle damping 70 to be the same between the supercon- 
ducting and the normal states, as was done in Onari et 
aVs previous study 3 -, they now take different values for 
the two states : 70 = 20meV for the normal state and 
7 = lOmeV for the superconducting state. The difference 
of the overall 7 form between the two papers is shown in 
Fig. [TJ The larger the difference of 7 between the super- 
conducting and normal states, the greater the enhance- 
ment (hump in the s ++ or the resonance peak in s±) of 
Im^s at q ~ (jr, 0) in the superconducting state over that 
in the normal state, in favor of their conclusion that the 
hump in the s ++ is large enough to explain the experi- 
ments while the peak in the s± is too large, as opposed 
to our conclusion. We do not deny the possibility that 
7 values may be different between the two states in the 
actual materials, but it is important to notice that the 
large difference among the studies lies here. 

On the other hand, Onari and Kontani claim that 
our conclusion 2 that "the hump in the s ++ is small for 
A = 5 ~ lOmeV compared to the case of A = 25meV or 
larger (as was taken in Onari et al's previous paper—)" 
is incorrect, based on the observation that Im% s in the 
region uj < 2A is finite, whereas it actually should be 
in a rigorous calculation. This fmiteness is in fact due 
to the approximation adopted in Onari et al.'s previous 
paper 3 . However, as shown in Fig. 4 of their new paper 1 . 
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FIG. 1: The difference in the quasiparticle damping between 
Onari et al. and Nagai et al.. 



the result of improving the uj < 2 A region only results 
in a suppression of the hump in the s++ state. There- 
fore, the improvement in the uj < 2A region does not 
affect our conclusion regarding the enhancement around 
uj = 2 A ~ 3 A in the s++. As for the s± state, we do 
admit that the broadness of the s± resonance peak in 
Fig. 1(a) of our paper may be due to the approximation, 
but actually we made no statement regarding this broad- 
ness, as opposed to what is mentioned in section IIIC of 
Onari-Kontani's paper. In fact, we showed another cal- 
culation in Fig. 1(b) of our paper— (which Onari-Kontani 
hardly mentions about), showing a sharper peak. As a 
common feature between the two cases, we mentioned 
only about the height of the peak in the s± state with 
respect to that in the normal state, which is highly de- 
pendent on the strength of the quasiparticle damping in 
the normal state. 

Also in Fig. 5(a) of refUJ, they refer to our calcula- 
tion result (Fig. 1(a) of our paper—) and mention that 
the peak position of the s ++ hump is close to uj = 2A, 
whereas it should be closer to 3 A. Here, they have failed 
to mention the value of U in our paper, which is actu- 
ally U = 1.375eV. The peak position is actually depen- 
dent on U, as already seen from the difference between 
U = 1.3eV and U — 1.32eV in Fig. 5(a) of their paper. 
We have performed a calculation for several values of U, 
and in fact the peak position systematically varies, as 
shown in Fig. [2] This variance on U could be to some 
extent due to the approximation adopted here, but in any 
case, the peak position is hardly relevant to our conclu- 
sion because in Fig. 1(b) of our calculation 2 , which Onari- 
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FIG. 2: ^/-dependence of Im^s at q = (7r,7r/16) in the s++- 
wave state with A = 5meV and y s = lOmeV. 
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FIG. 3: Imxs at q = (tt, tt/16) with [/ = 1.3eV. Blue triangle 
is for the s++ state with A = 5meV and 7 S =10meV. Red 
circles and green squares are for the normal state with the 
damping parameter 70 = lOmeV and 20meV, respectively. 



Kontani makes no comment on, the peak position is in 
fact around oj — 3 A. 

To further strengthen this point, we once again show 
in the present Comment, now adopting U = 1.3eV and 
A = 5meV, that our conclusion holds as far as we adopt 
the same values of the quasiparticle damping between the 
normal and the superconducting states. In Fig.[3j we see 
that the enhancement in the s++ state above the nomral 
state is small when we adopt 7 = lOmeV for both states. 
On the other hand, if we adopt 7 = 20meV only for the 
normal state, the s++ hump now appears to be large, con- 
sistent with Onari-Kontani's claim. However, there does 
lie a quantitative difference between their calculation and 
ours : in the high energy regime of u = 6A or more, Imx s 
in the superconducting state with 7 S = lOmeV merges 
to the normal state values with a different 7o(=20meV) 
in Onari-Kontani's calculation^ (Fig. 2 and Fig. 3 of their 
paper), while in our calculation it merges to the normal 
state values with the same 7o(=10meV). In the high en- 
ergy regime, the effect of the superconducting gap should 
disappear, so that we believe our results are more physi- 
cal. We are not sure about the origin of this discrepancy. 

To conclude, the main difference between Onari- 
Kontani's calculation^ and ours (especially Fig.l(b))^ lies 
in the choice of the quasiparticle damping in the normal 
state, not in the method or the accuracy of the calcu- 
lation. Their claim that the conclusions of our paper 
are "incorrect based on inaccurate numerical calculation" 
has no basis and is in fact not correct. If we consider the 
possibility of the 7 values being different between the su- 
perconducting and the normal states (which we do not 
deny), then the ambiguity of determining the supercon- 
ducting gap form solely from the neutron scattering at 
~ (71", 0) will be greater. Then our proposal to look into 
wave vectors other than (n, 0) will have even more im- 
portant implication^. 
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